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In this note the ILC's capabilities for detecting WIMPs and measuring their properties 
are studied. The expected signal cross section is derived in a model-independent way 
' from the observed relic density of Dark Matter. Signal events are detected by means of 

. initial state radiation (ISR). The study is performed with a full simulation of the ILD 

detector. The results show that WIMPs are observable at the ILC if their coupling 
to electrons is not too small (O(0. 1)). Their masses can be measured with a precision 
of 1 to 2 GeV. The accessible phase space can be increased significantly using polarised 
beams, especially if the positrons are polarised as well. 



1 Introduction 



Due to its clean experimental environment the International Linear Collider (ILC) offers the 
possibility to look for Weakly Interacting Massive Particles (WIMPs) in a model-independent 
■ way. This has been pointed out in ^ [3] . In this study the expected sensitivity to such a 

Q\ \ WIMP signal, the achievable mass resolution and the influence of beam polarisation are 

investigated using a full detector simulation of the International Large Detector (ILD) of 
^" ■ the ILC. 



2 Model-independent WIMP production cross section 



It is assumed that the cosmic Dark Matter (DM) component is only due to one new type 
of particle, and that the relic density is determined by pair annihilation of WIMPs into 
5^ ■ SM particles. These annihilations may produce an e + e~ pair (xx ~ > e + e _ ) with an un- 

known branching fraction K e , where K e > 0.3 is strongly motivated by recent PAMELA 
results [U [5] . The total annihilation cross section <r an is then determined by the observed 
DM density Qdm- The annihilation cross section required to match the observed relic den- 
sity is in the order of a few pb [3] . Using crossing relations one can derive the expected cross 
section for the reverse process, i. e. e + e~ — ► xXi which could be observable at the ILC. The 
resulting cross section contains as free parameters: 

• the e + e~ branching fraction K e 

• the mass of the WIMP M x and its spin S x 



the angular momentum Jo of the dominant partial wave (Jo corresponds to the spin 
of the exchange particle in the annihilation (production) process). 
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Since the produced WIMPs leave the detector without any further interaction an additional 
photon from ISR is required (e + e~ — > xXl)- The main background is the radiative neu- 
trino pair production e + e~ — ► vvy. This reaction is mediated dominantly by t-channcl 
VF-exchange if the center of mass energy is significantly above the Z°-pole. Therefore 
the vvj background is strongly reducible depending on the polarisation state of the ini- 
tial electrons and positrons. 

In this model-independent approach no constraints are imposed on the hclicity structure of 
the WIMP couplings to electrons, which therefore enter the production cross section as a 
new free parameter. We consider three different scenarios for the helicity structure of the 
WIMP couplings: 

• the same as the SM charged current interaction, i. e. only fc(e^e^) is nonzero 

• parity and helicity are conserved: ft(e^e^) = K(e^eJ) 

• opposite to the SM charged current interaction, only ^(e^ej) is nonzero. 

In the last two cases a significant enhancement of the signal to background ratio is expected 
if electrons (and positrons) beams are appropriatly polarised. 

3 Software and reconstruction tools 

For the SM vvy background the SLAC SM Whizard event sample [6] is used. This sample 
is generated with the Whizard event generator including beamstrahlung and beam energy 
spread from Guinea Pig for nominal ILC beam parameters. Each e + e~ — > wy event con- 
tains up to two additional ISR photons. The sample size corresponds to 10 fb _1 of luminosity 
at y/s = 500 GeV. The background sample is reweighted with respect to the energy and 
polar angle of the detected photon according to the WIMP production and background 
cross sections. The benefit of this method is that a full signal sample can be obtained for all 
investigated combinations of cross section parameters with only one simulation and recon- 
struction cycle, reducing the amount of processing time significantly. Since the predicted 
signal cross section assumes only one photon to be radiated, the center of mass energy of 
the hard subprocess is used for the application of weights. 

The event sample is then fed into the full Mokka 06-06-p03 ILD detector simulation (de- 
tector model LDCPrime_02Sc) with a 3.5 Tesla magnetic field, and reconstructed with 
MarlinReco using the Pandora Particle Flow algorithm (PFlow) [7]. 

The event simulation and reconstruction process introduces a smearing of the recon- 
structed photon candidate energies due to the intrinsic energy resolution of the detector and 
the performance of the Particle Flow algorithm. Figure Q^a) shows the energy distribution 
of the most energetic photon candidate identified by the Pandora algorithm. The radiative 
return to the Z° at ~ 240 GeV is with a width of about 6.4 GeV significantly broadened 
compared to the spectrum on generator level (Fig. HJb)). 

4 Cluster splitting 

The full reconstruction yields less photons at high energies resulting in a shift of the Z° 
return from the expected central value of 241 GeV to 237.7 GeV (Fig. [Ha)). This behaviour 
is quantified in Figure[2]Ja) showing the average ratio E rec / E gen of the reconstructed photon 
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Figure 1: Energy spectrum of the most energetic photon of the i/wy background (a) after 
reconstruction using Pandora PFlow and (b) on generator level. 
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Figure 2: (a) Average ratio of reconstructed photon energy to generated photon energy 
before and after merging, (b) Energy spectra before and after merging procedure. 



candidate's energy to the photon energy on generator level. The Pandora PFlow algorithm 
tends to split photon clusters stemming from one generated photon into several photon 
candidates. This lowers the ratio E rec /E gen especially at high energies (black triangles). 
This effect can be compensated by applying a simple merging procedure. Neighbouring 
photon candidates are combined to form a new set of photon candidates, with which an 
average ratio of E rec / E gen sw 1 is recovered (red squares) . The Z° return is also shifted back 
to its expected value (see Figure HJb)). 



5 Preliminary analysis results 

So far the following scenarios have been investigated: 

• WIMP spin: p-wave annihilation (Jo = 1) for S x = 1 and S x = 1/2 

• WIMP couplings: K(e^e^) > 0,«:(e^ej) > and Ac(e^e^) = re(e^e£) > 

• Polarisation: unpolarised beams, e~ polarisation only (P e - = 0.8) and additional e 
polarisation (P e - = 0.8 and P e + = 0.6). 
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Figure 3: Observational reach (3tr) of the ILC for a Spin-1 WIMP in terms of WIMP mass 
and K e for three different chirality and couling scenarios. Full line: P e - = P e + = 0, dotted 
line: P e - = 0.8, P e + = 0, dashed line : P e - = 0.8, P e+ = 0.6. 



5.1 Observational reach 



For each combination of these parameters the ILC reach for a 3er observation has been deter- 
mined as a function of the WIMP mass M x . The reach has been evaluated for an integrated 
luminosity of 500 fb _1 at -y/s = 500 GeV using a larger event sample from a previous de- 
tector simulation [2]. Due to the high SM neutrino production background, the sensitivity 
has been obtained statistically by using fractional event counting [8] as implemented in the 
ROOT class TLimit. Figure |3] shows the expected ILC sensitivity for Spin-1 WIMPs in terms 
of the minimal observable branching fraction to electrons function of the WIMP 

mass. The regions above the curves are accessible, where the full line gives the result for 
unpolarised beams, the dotted line for P e - = 0.8 and the dashed line for P e - = 0.8 and 
P e + = 0.6. Figure[3^a) shows the case where the WIMPs couple only to lefthanded electrons 
and righthanded positrons, fc(e^e^), Fig[3Jb) shows the parity and helicity conserving case, 
K ( e Z e fl) = K ( e R e t)> while Fig [3Jc) right plot is dedicated to the case that the WIMPs 
couple to righthanded electrons and lefthanded positrons (^(e^ej)). 

For the latter two scenarios polarised 
beams increase the reach significantly, es- 
pecially the additional positron polarisation 
increases the accessible range in K e by about 
a factor of two. In case of a Spin-i WIMP 
the sensitivity is somewhat worse, but again 
beam polarisation extends the observable 
part of the parameter space significantly pQ . 
In all presented cases the ILC is sensitive to 
the branching fraction K e over a large range 
of possible WIMP masses and down to val- 
ues below 0.3 indicated by the PAMELA 
results. 
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Figure 4: Recoil mass 
180 GeV Spin-1 WIMP. 
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5.2 Mass resolution 

If WIMPs are observed at the ILC, their mass can be determined from the recoil mass 
distribution of the photons: 

^recoil = S - 2VIE 7 (1) 

Figured] shows an example for the recoil mass distribution for a 180 GeV Spin-1 WIMP with 
both beams polarised. The WIMP signal shown in red (dark grey) on top of the SM neutrino 
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Figure 5: A% 2 from mass determination by a template method for a Spin-1 WIMP with 
A/x=150 GeV. (a) parity and helicity conserving couplings, (b) K(e^eJ). Full line: P e - = 
P e + = 0, dotted line: P e - = 0.8, P e + = 0, dashed line : P e - = 0.8, P e + = 0.6. 

background kicks in at M reC oi] = 360 GeV. From this distribution the WIMP mass can be 
reconstructed e. g. with a template method. For this procedure, only 200 fb _1 of the available 
MC sample have been analysed as dataset, the rest is used for the templates. Figure [3] shows 
the obtained A% 2 as function of the reconstructed WIMP mass for a 150 GeV Spin-1 WIMP 
for K e = 0.3. Fig.[5ja) shows the helicity and parity conserving case, while (b) illustrates the 
case that the WIMPs couple to righthanded electrons and lefthanded positrons, ft(e^ej). 
Again the full line gives the result for unpolarised beams, the dotted line for P e - = 0.8 and 
the dashed line for P e - = 0.8 and P e + = 0.6. Without any beam polarisation, the mass 
resolution is about 4 GeV, which is reduced to about 1.2 GeV by switching on the electron 
polarisation. Additional positron polarisation improves the resolution by another factor of 
two to about 0.6 GeV. 

Similar results are obtained for Spin-i WIMPs pp. As for the observation reach, the 
situation is slightly worse than in the Spin-1 case, but again the use of beam polarisation 
leads to a significant gain in resolution. 

6 Conclusions and outlook 

This study is one example of a physics analysis using the full simulation of the ILD detector 
concept. It is part of the optimisation effort for the ILD detector, and is further improved 
in parallel to the ongoing optimisation of the detector model. It demonstrates that there is 
a good chance of detecting WIMPs at the ILC in a model-independent way, and to measure 
their mass with a precision of about 1 GeV. The preliminary results show that the reach of 
the ILC on the branching fraction K e of WIMP pair annihilation to electrons covers values 
of n e down to ~ 0.01, which is well below the n e > 0.3 indicated by the recent PAMELA 
data. Both the range in phase space, as well as the mass resolution improve significantly 
when polarsised beams are assumed. Typically the use of 80% electron polarisation gives 
improvements of a factor of two over unpolarised beams, whereas an additional positron 
polarisation of 60% yields another factor of two. 

The results will be updated using different Particle Flow algorithms and Photon finders. 
Reducible backgrounds as well as beamstrahlung will be included. The obtained results will 
also be compared to specific SUSY scenarios in which the only open SUSY channel at the 
ILC's center of mass energy is radiative neutralino production [9]. 
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